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^ ■ Abstract 
m 

^ . We point out that the correlation between folding times and a = {Tq 

o 
o 

. and folding transition temperatures was already established in 1993 before 

0^ ■ the other presumed equivalent criterion (folding times correlating with Tj 



Tf)/TQ in protein- like heteropolymer models where Tg and Tf are the collapse 



alone) was suggested. We argue that the folding times for these models show 



. no useful correlation with the energy gap even if restricted to the ensemble of 

c : 

O ■ compact structures as suggested by Karplus and Shakhnovich. 

O 

> 

•i-H 

' In a recent article ^ we (KT) showed that for certain lattice models of proteins the 

folding times, rj, correlate extremely well with a = (Tq — Tf)/Tg where Tg and Tf are well 
defined thermodynamic collapse and folding transition temperatures, respectively. We also 
demonstrated that there is "no useful correlation between folding times and the energy gap 
between the native conformation and the first excited state" . In response to these results 
Karplus and Shakhnovich (KS) try to argue (i) that the folding criterion used by KT is 
"essentially the same as one introduced earlier 0,^" and (ii) that the energy gap used by 
KT is not " appropriate" . We take up these two issues separately. In addition, we also show 
that folding times do not correlate with the energy gap Acs restricted to the ensemble of 
compact structures as KS desire. 
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First we settle the historical claims to priority. The criterion used by KT was already 
established in 1993 in which was published before 0] was even received for review. There 
is no plot of folding time versus T/ in 0, an article that was submitted after P was accepted 
for publication. It was shown in ^ (cited in |^) that folding times correlate with a. We 
showed using numerical results and theoretical arguments that "there appears to be useful 
correlation between folding time r,. and a = l — Tj/Tg: the smaller the value of a the smaller 
the value of r^" p. 

Having addressed the historical claim to priority we now examine the statement that our 
folding criterion is "essentially the same as one introduced earlier" This claim is based 
on the observation that there is a correlation between folding time and Tf seen in "right 
most part of Fig. (8a) of [|]". Fig. (8a) of 0| shows the folding time for 10 sequences of 
which 5 lie in the "right most part". The simulation temperature for the data in Fig. (8a) 
is Tg = 1.0 and the Tf for all the sequences lie in the range 0.63 < Tf < 1.06 (see Table 
I in Q). Of the 5 sequences that lie in the "right most part" of Fig. (8a) only one has 
Tf > Tg. The apparent correlation (for the five sequences for which the folding times merely 
change by a factor of 3) claimed by KS occurs when the native state is not stable. In Fig. 
(8b) in 1^ a plot of folding time as a function of at = 0.7 (which is below Tf for all 
sequences except two) is shown. In this figure, with Tf > Ts, one observes no correlation 
whatsoever between Tf and Tf. In contrast, our results [jll show that Tf (which spans six 
orders of magnitude) correlates well with a. It is the interplay between the two intrinsic 
sequence dependent temperatures Tg and Tf that seems to correlate with Tf under conditions 
when the native state has the highest occupation probability. 

It was forcefully stated in that the "necessary and sufficient" condition for folding 
in these models is that there should be "a pronounced energy gap between the native and 
first excited state for the fully compact ensemble" The reason we did not display Tf 
as a function of Acs in is that roughly half of the sequences had non-compact native 
conformations. It is, therefore, not clear why one should be restricted only to the ensemble 
of compact structures. However, the folding times for sequences in [|] when plotted as a 
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function of Acs, as KS desire, also does not show any correlation (see Fig. (la)). There 
are a large number of sequences, all with roughly the same folding times, but with very 
different values of Acs- The lack of correlation between the energy gap and foldability has 
been established for other models as well P-[Tl[ . 



There are a few additional comments that are worth making. (1) KS claim that Tf should 
correlate well with both Tf and Acs. This implies that Tf should correlate with Acs. In 
Fig. (lb) we plot Tf as a function of Acs- The lack of correlation is evident. (2) Unlike Tq 
and Tf, Acs is not experimentally measurable. Moreover, with the exception of two studies 
|[TO| , P^ , the energy gap is not easily computable for simple off-lattice models. Thus the 
practical utility of Acs is not clear. (3) We should emphasize that all the criteria proposed 
for folding times should be viewed as statistical. This means that, under conditions when 
the native state is significantly populated, Tf may correlate with certain properties intrinsic 
to the sequence provided a number of sequences is studied. Our results suggest that a could 
be one of such properties. 
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FIGURES 

Fig. 1. (a) Folding time t/ versus Acs- (b) as a function of Acs- Solid circles 
correspond to 15-mer and open circles are for 27-mer. Right and upper axes are for 27-mer. 
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